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TECHNICAL NOTE 3458 

UNSTABLE CONVECTION IN VERTICAL CHANNELS WITH HEATIHG 
FROM BELOW, INCLUDING EFPECHS OF HEAT 
SOURCES AND FRICTIONAL HEATING 
By Simon Ostrach. 


SUMMARY 

The effects of heat soiorces and frictional heating on the laminar 
fully developed channel flow subject to a body force between two paral- 
lel plates oriented in the direction of the body force are analyzed. 
Solutions are obtained for combined forced- and natural- convection flows 
for the cases in which the wall temperature variations are linear and 
(l) the wall teiaperatures are specified, (2) the walls are both insu- 
lated, and ( 3 ) the net mass flow in the cbs.nnal is zero. These solu- 
tions depend on the Rayleigh number which was previously fovmd to be 
the factor determi nin g the stability and type of flow for horizontal 
and vertical layers of fluids heated from below but without heat sotrrces • 
or frictional heating. Similar stability characteristics are displayed 
in the present problem, and the heat sources affect the flows only in a 
quantitative manner. 

When the effects of frictional heating are considered, two distinct 
solutions are obtained for each set of parametric values. Solutions 
neglecting frictional heating correspond to a different exact solution 
for values of the Rayleigh number smaller and larger than the critical. 
The related approximate and exact solutions are essentially coincident 
when the frictional-heating parameter is rthbI l but differ for unit order 
values of this parameter. The approximate solutions (l.e., those neg- 
lecting frictional heating) are shown to be always invalid for Raylei^ 
numbers near critical. 


INTRODUCTION 

CJonsiderable theoretical work (refs. 1 to 5) has recently been done 
on problems dealing with internal natural-convection flows because of 
the many new practical applications of this phenomenon. The heat- 
transfer results of reference 5 are experimentally verified In reference 
6 . All this work is concerned with flows in vertical enclosxires . 
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One of the many interesting aspects of the natural-convection proc- 
ess^ which was originally studied theoretically by Eaylel^ in 1916 
(ref. 7) for horizontal rather than vertical layers of fluid, is the in- 
stability of the flow associated with heating a fluid at rest from be- 
low. When a stationary fluid iias some layer with a density greater than 
others lying below it, its equillbritim becomes iinstable in the sense that 
even a small disturbance may result in a completely altered regime. 

Benard (ref. 8} and others had previously shown estperimentally that such 
a configuration leads to a cellular type of flow. Numerous theoretical 
and experimental investigations have been performed to substantiate, 
clarify, and exctend the early findings, but these works were restricted 
to horizontal layers of fliild. lEhese studies, which are briefly summa- 
rized in reference 9, were essentially of the stability type in which 
eigen relations of the physical and cell- shape parameters were deter- 
mined for neutral stability, and the details of the flow were usually 
not investigated further. 

The results of this work on horizontal layers indicate that the 
flow regimes depend on a dimensionless parameter sometimes referred 
to as the Rayleigh number, vdiich is the product of the Prandtl and 
Grashof mmibers. It has been demonstrated experimentally that there 
are essentially two different laminar types of flow: (l) a columnar 

type (ref. 10), occurring after a relatively low value of the Rayleigh 
number is attained, and (2) the cellular type Yobserved first), occur- 
ring for larger values of the Rayleigh number -'^approximately 1700 for 
two horizontal parallel plates). . 

— In reference 5, it was pointed out that heating from below might be 
expected,jb©-affect the stability of vertical as well as horizontal lay- 
eja-dSspite the differences between the two configuratlonB . As a matter 
of fact, heating from below In vertical enclosures is often actually en- 
countered; a discussion of such a configuration in an atomic power vinlt 
is presented in reference 11. Configurations of this type can also be 
found in mai^ other fields besides that of atomic power; for exanple, 
this phenomenon can also be found in conjunction with the cooling of 
turbine blades by natural convection. (The writer is unaware of any 
previous presentation of this connection. ) Accordingly, the problem 
treated in reference 5 (namely, the fully developed flow with temper- 
atiore increasing linearly upward in a vertical channel) was reconsid- 
ered in reference 12, with the modification that the heating was from 
below; that is, the temperature was specified to decrease linearly 
upward. 

In reference 12, it was found that the stability characteristics 
and the critical values of the Raylei^ number for the vertical fluid 
layers were similar to those for horizontal layers. Representative 
velocity and temperature distributions for various sets of boundary 
conditions were determined from explicit expressions in reference 12 
neglecting frictional heating. It was shown that, in certain ranges 
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of the parametric values, the flow and heat transfer associated with 
heating frcan helow are appreciably different trcm the corresponding 
quantities for the identical configuration in idiich heating is not 
from below. 

The analysis of reference 12 is extended herein to include the ef- 
fects of uniform internal heat sources in the fluid and of frictional 
heating. The iiaportance of frictional heating in the natural- convection 
phenomenon was first discussed in reference 1 and later in references 5 
and 13. GThe solutions obtained with frictional heating are exact, so 
that the assxm^tion of the neg31gibility of frictional heating in refer- 
ence 12 can be checked. 


AMLYSIS 
Basic Equations 

Ihe fully developed laminar flow of visoous fluids subject to a 
body force between two plane vertical parallel surfaces open at both 
ends (see fig. l) is considered herein, as in references 1 wnd 5 . 

(More generally, the surfaces could be taken to be oriented in the di- 
rection parallel to any generating body force . ) In contradistinction 
to the work in references 1 and 5, the fluid here is heated from below, 
so that the axial (vertical) temperature gradient is specified as nega- 
tive, It is further assumed that the physical properties of the fluids 
are constants, except that the essential influences of density changes 
on the flow are taken into account insofar as they modify the effects 
of the body forces . This last assumption is usually employed in prob- 
lems of this nature , . and its Justification is discussed in references 
1 and 14. 

Ikider the conditions stated (which are essentially those for 
Poiseuille flows together with the specific thermal conditions), the 
temperatin-e can be expressed as the sum of a linear function of the 
vertical coordinate and arbitrary function of the horizontal coordinate: 

T* = AX + T(Y) ( 1 ) 

Further, the continuity equation is identically satisfied, and the mo- 
mentum and energy equations in dimensionless focrm are 

u'‘ + T = -CK ( 2 ) 

t" + Ra u + (u')^ + ctJC = 0 (3) 

(All symbols are defined in the appendix.) The primes denote differen- 
tiation with respect to the dimensionless horizontal coordinate y. The 
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terms in equation (2)(fi‘om left to right, respectively) denote the vis- 
cous, huosrancy, and axial-pressure forces, and the terms in equation (s) 
denote the conduction, convection, frictiooal heating, and heat sottrces. 
!Btie details of. the development of these equations are essentially pre- 
sented in reference 5, the only modification in the present work being 
that the longitudinal (vertical) teriperature gradient =* A is 

taken as negative- Note that the Rayleigh and Grashof numbers are modi- 
fied herein in that the characteristic temperature is Ad rather than a 
temperat\3re difference as is conventional. 


The constant C ■which appears as a parameter in the problem merely 
specifies the tenperatxjre level (see eq. ^2)) and is defined by (see 
ref. 5) 



This parameter must in some way be related to the physics of the prob- 
lem. Rrom equation (4) it can be seen that C could be determined from 
the longitudinal pressure gradient; -that is, C is essentially connected 
■with the end conditions to which the channel is subject. The parameter 
C can also be related to the end conditions by the mass flow in the 
channel, which rexoains Invariant over the entire channel length- Such 
a relation will be developed sxibsequently. The solutions of equations 
(2) and (3), together ■with appropriate boundary conditions, will, 
through the parameter C, apply to both natural-convection and ccmbined 
natural- and forced-convection flows. The forced-convection pressxnre 
gradient merely alters the magnitude of C- 

The system of equations (2) and (3) can, by eliminating T, be 
■written as a single fourth-order nonlinear ordinary differential 
equation: 

u^"'^ - Ra u - (u')^ - (xK => 0 (s) 

Equations of the same form as in reference 5 are ob-tained by letting 


V 


64 


(Ra u + oK) 

(6) 

2y - 1 

(7) 

II 

(8) 




V 


i< 
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where equation (?) is merely for convenience. Equation (s) 500011108, for 


viv _ = 0 


( 9 ) 


OJ 

U1 

-J 

w 




Boundary Conditions 


The present problem will be considered for two different sets of 
theimal-boundary conditions at the channel surfaces, namely, (l) that 
the wall teniperatures are specified and (2) that both walls are insu- 
lated. Another case, which was essentially a combination of the above 
conditions, was treated in reference 12 } but, since no basically new 
results were indicated, this case is not discussed herein. 

For both sets of thermal conditions, the no- slip condition of vis- 
cous fluids must be satisfied; that is, 

v(-l) = v(l) = A (10) 


where 


A 


0^ 

64 


( 11 ) 


Specified side-wall temperatures . - For the first set of thermal 
conditions the side-wall temperatures are specified. The ten5>erature 
varies linearly along the walls (equal slopes on both walls), but the 
walls may each be at different teanperatirres . Thus, equation fl) im- 
plies uniform heat fl-ux across each surface. When equations ( 2 ;, (s), 
and (7) are used, this condition can be written as 


v"(-l) = 
v"(l) = nJ's/R 


where 


J = 


RaV^CK 

64 


( 12 ) 

(13) 

( 14 ) 


and 


''w-i 


CC-A)d 


n = 1 + 


(15) 
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■where n Is essentially a meastire of the difference hettreen the two 
wall temperatures. 


Both side •walls ins'ulated . - In order to simulate more closely the 
conditions used In the study of the horizontal layers where the side 
walls were taken to be insulated surfaces of symmetry in the fluids the 
second set of thennal conditions chosen herein is that both side walls 
are Insulated; that is^ 


or 





v'"(-l) = v'"(l) = 0 


(16) 


Solutions 

The nonlinear term in equation (9) is associated with frictional 
heating. If the effects of frictional heating are neglected^ the prob- 
lem then consists in solving the equation 

- Rvq = 0 (17) 

together -with the proper boimdary conditions. For no internal heat 
sources (a = o) , this problem is solved in reference 12 for three 
sets of boundary conditions. The solutions of equation (l7), identi- 
fied by the subscript zero, are hereinafter referred to as the approx- 
imate solutions. The effects of a uniform distribution of heat sources 
(a ^ O) are determined herein. 

In the present paper the effects of frictional heating are deter- 
mined from a direct solution of the complete equation (s), by means of 
a forward integration technique on a high-speed computing machine (Card- 
Programmed Calciilator) 

The effects of frictional heating are clearly discernible by com- 
paring the solutions of the complete and the approximate equations (eqs . 
(5) and ( 17 ), respectively). The solution of eqiiation (l7) also repre- 
sents an extension (including heat sovrrces) of the work in reference 12. 
The general solution of equation (l7) is of the form 

vq = Cl cos + Cg sin riA,, + C3 cosh E^/'^ti + C4 slnh E^Ati 


4he effects of frictional heating co\ild also be determined by an 
iteration procedure as described in reference 5, however, this procedure 
■will not be used herein. 
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For each given set of boimdary conditions^ the solutions are presented 
in the succeeding sections in terms of Vq, and the corresponding tem- 
perature differences t can he obtained by use of equation (2). 


Specified side--wa.1 1 temperatures - The solution of equation (l7) 
subject to the boundary conditions (eqs. (lO)^ (l2), and (l3)) is 


. . J(n - l) . Jfn + 1) 

■'^0 ^ ~ 2 ^ "'^01 + — 2 ^ ^02 

where 

1/cosh , cos 

2Vcosh rIA cos rIA J 

l/sinh E^Ati sin rV\ J 

"" 2\sinh rIA ~ sin rV^ / 

_ l /cosh Rl/S . COS 

2^cosh RIA cob r1/ ^ / 


(18) 


(19) 

( 20 ) 
( 21 ) 


The solution given by equation (18) is a function of the parameter 
C. In order to relate this solution to a given physical problem, C can 
either be determined from the channel presstire gradient through equation 
( 4 ) or from the dimensionless mass flow through the eacpression 


M = UQ dy (22) 

Equation (22) is particularly convenient in the case of no net mass flow 
through the channel, for then C can be directly related to the param- 
eters a, n, and Ea. 

Examination of the solution neglecting frictional heating (eqs. 

(I 8 ) to ( 21 )) shows that there are critical values of the Rayleigh 
number ^ for which the solutions become infinite. These values, 
namely, Rajj- = (lor)-, where Is: denotes integers, obtained herein con- 

sidering heat som'ces are identical with those given in reference 12 
for no heat soxjrces. It is significant that these critical values of 
the Eaylei^ number when k = 1 and 2 (Ra 2 _ = 97.41 and Ea 2 = 1558.55) 
correlate closely with the critical values found for horizontal layers 
heated from below (see ref. 9). Further discussion of all the solutions 
and their physical slmpll cations is given in the next section. 
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Both side walls Insialated . - The solution of the houniaiy- value 
prohlem specified by equations ill), (lO) ^ and (is) is 


^ _ xfsin R^/^cosh - sihh B^/^cos 

^ sin RV^cosh rV^ - sinh R^A cos r1/4 


(23) 


The critical values of Ra for equation (23) are those which satisfy 
tan R^A = tanh R^A or Ra^^^ S'jj^k + k 1, Ra^_ = 3803.22. 

Note that the boundary-value problem as s^tated for this case (both 
walls insulated) by equations (s), (s), (lO), and (l6) defines t to 
within an arbitrary constant. Therefore^ in order that t vanish at 
the wall (y = O) , which it must by its definition^ C must be so speci- 
fied as to accomplish this. Hence, for this case, t is determined 
from (see eq. (2)) t = u"(o) - u" . 

Solutions for special case simulating an enclosed channel . - In or- 
der to simulate, a completely enclosed channel, that is, one in which the 
ends are closed, the net mass flow as given by equation (22) can be 
specified as zero. For the case of linear side-wall temperatirres, no 
net mass flow will be obtained if 


C(n + l) _ a tanh rV^ + tan rV^ - 2rV^ 
2 4 rV 2 tanh R^/^- - tan R^/^ 


(24) 


For the case wherein both walls are Insulated, the condition of no 
net mass flow can be obtained only if there is no heat due to heat 
sources (a = O) . The sittiation for a =* 0 is as described in reference 
12; that is, no flow except for certain specific values of Ra (those 


satisfying 
ref. 12). 


cos Ra 


^/^cosh Ra^/^ =s 1 


or Rai 


■k 




as given in 


I 


V 


49 


MECHANICAL ANALOGY 

Equation (l7), obtained from the mathematical formulation of the 
problem treated herein, is identical with that describing the vibra- 
tions of uniform beams or of rotating shafts (see ref. 15). The case 
wherein both walls are insulated is analogous to the vibration problem 
if the ends of the shaft are fixed with no shear. The critical values 
of the Rayleigh nvunber correspond to the critical or whirling speeds 
of the shaft. 
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If the shaft has initial deflections or "bending moments at the 
endS; the rotating shaft "becomes dynamically unstable and the solu- 
tions become infinite at the critical speeds. This situation has its 
counterpart in the present problem 'wherein the side-'wall temperatures 
are specified. 


EESULTS AMD DISCUSSION 


The discussion of the results is divided into t'wo main sections 
according to the thermal "boundary conditions. The main sections are 
further subdi-vided by first neglecting the frictional heating effects 
and then including them in order to demonstrate clearly the effects 
of heat sources and of frictional heating. 

The primary effect of heating from below, namely, the instability 
of the flow, is demonstrated by comparing the present resiilts with those 
for the sta"ble case (A > O) reported in reference 5. In -this regard, it 
should be noted that, since CK, oK, CA, and K/A are all independent 
of A, the problem treated herein (for A < O) differs frcm that in ref- 
erence 5 (for A > O) only in the sign of -fche convection term. There- 
fore, in order to canpare the two cases properly, the boundary condi- 
tions must be identical in the two problems. Choosing C wi'th opposite 
signs in the two problems makes them identical. Since C = -1 for the 
present calculations, the corresponding cases for A > 0 are computed 
■with the use of reference 5 for C =« 1. 


The resiilts are all presented in dimensionless form, -v/ith the phys- 
ical quantities related to the dimensionless ones by 



T = yd 


(27) 


Linear Side-Wall Temperatures 

Mo frictional heating . - In order to investigate closely the flow 
and heat transfer, including the effects of heat soiorces for values of 
Ra different from 'the critical, calculations of -the velocity ar>d tem- 
perature dl'stributions were made from the approximate solutions for 
a = 10, K = 10, and the t'wo side -walls ei'ther at the same temperature 
(n » l) or at different temperatures (n = 0, for convenience) and 
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Ra = 10^ 81^ 100 j 1000 j and 1600, so that the critical values of the 
Rayleigh number are straddled. The calculations of reference 12 for no 
heat sources (a = O) were also extended herein over the larger range of 
Rayleigh number values . All these ccmputations are presented in figures 
2 and 3 for the symmetric case (n =» l) and in figures 4 and 5 for the 
asymmetric case (n = O) end are labeled with Uq and Tq to denote 
that they are the solutions obtained neglecting frictional heating. 

From the solutions (eq.s. (iS) to (2l)), it can be seen that the 
only critical value of ^ for n = 1 in the range 0 < Ra 4: 1600 
(for which the calculations were made) is Ra = jr 4 =* 97.41. Com- 
parison of the velocity profiles (fig. 2) for each fixed a and both 
walls at the same temperature shows a different flow direction on ei- 
ther side of this critical point, the distributions remaining essen- 
tially the same shape except for boundary- layer effects as Ra is 
further increased to 1600. This behavior can be clearly observed in 
figure 6, -vdiere the velocity extrema are presented as functions of 
Ra for a ?s 0. The temperature distributions (see fig. 3) for a = 0 
are of opposite cinrvattire on either side of Ra = 97.41; but for a = 10 
there is no such change, but merely a gradual decrease of the maximum 
values followed by an increase for large Rayleigh numbers. 

For the walls at different temperatures (n = 0), the two critical 
values of Ra in the range up to 1600 are * 97.41 and 
(2rc)^ =» 1558.55. For this case, it can be seen from figures 4(a) and 
(b) that again the velocity changes direction across Ra =* 97.41 for 
a = 0 and a = 10. For values of Ra > 100, the flow is no longer 
entirely in one direction but is in opposite directions in adjacent 
sections of the channel (see fig. 4(c)). Beyond the second critical 
point, Ra = (2x)^ = 1558.55, the flow directions in the adjacent sec- 
tions are reversed. The temperature distributions (fig. 5) are also 
appreciably altered in passing the critical values of the Rayleigh 
number . 

The effects of heat sources in the fluid can be determined by com- 
paring the curves for a = 0 with those for a = 10 in figures 2 to 5. 
For Ra = 10 (figs. 2(a) and 4(a)), the addition of the heat (a = lo) 
for these computations) causes the flow direction to be reversed essen- 
tially (with different magnitudes). Near the first critical point, the 
heat sources greatly reduce the velocity magnitudes . The addition of 
heat sources, of course, alters the temperature distributions and, 
hence, alters the heat transfer greatly. 

In order to study the effect of heating from below, calculations 
made for the stable case (A > O) frcrn reference 5 are superposed on 
figures 2 to 5. The boundary conditions are specified to be identical 
for both problems, so that the two cases correspond except for the 
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heating from, helow. Ccai 5 )arison of the velocity and temperatiare profiles 
shows that the heating from below changes the flow and heat-transfer di- 
rections- and magnitudes in seme ranges of Ra values. In all cases j 
greater velocities result -with heating from, below. Since ^ as has been 
previously discussed, the heat so-urces can also change the flow and 
heat -transfer magnitudes and directions frem the a = 0 case, the ef- 
fect of the heat sources in this ccoparison is merely to shift seme of 
the ranges of Ra values for which the stable and unstable flows and 
heat transfers are unidirectional. 

Representative velocity and ten^ierature profiles for the case simu- 
lating a fully enclosed region (i.e., zero net mass flow) are presented 
in figures 7 and 8 for Raylei^ numbers of 10, 100, and 1600 and for 
heat-so-urce parameters of 0 and 10. For small Ra (see figs. 7(a) and 
8(a)), the -velocity profile with no heat sources (a = o) is antisymmet- 
ric and the tenrpera-ture profile is essentially linear. For heat sources 
in the flioid (a = lO) , the -velocity profile is symmetric about the chan- 
nel axis and the temperature profile shows the effects of the Increased 
convection by its almost parabolic shape. For large Rayleigh numbers 
(figs. 7(c) and 8(c)), both the velocity and temperature profiles for 
both a = 0 and a = 10 are antisymmetric, but for each different a 
the directions are opposite. Further, since Ra = 1558.55 is the first 
critical value due to heating from below foi* this case, the direction 
change fron that for small Ra is also clearly e-vident. Note that the 
parameter n for no 3iet mass flow and fixed C changes with a, and 
the appropriate values are indicated in the figures. 

Frictional heating . - In reference 5, it -was poin-bed out that, in 
general, the effects of frictional heating could be important if -the ra- 
tio K/Ra is of -unit order or larger. Iherefore, in order to see how 
well the solutions neglecting frictional heating approxlma-te the exact 
solutions (eq. (5)), the conplete equation including the effects of 
frictional heating was solved for several cases by means of a Card- 
Programmed Calculator. A forward integration technique similar to that 
described in appendix B of reference 16 was employed. The results of 
these calculations, denoted by u and t -with no subscripts but -vd.th 
superscripts to indicate the exact solutions, are presented in tables I 
and II and also in figures 2(a) and (d), 3(a) and (d), 4(a) to (d), and 
5(a) to (d) . The significance of the superscripts -v/ill be discussed 
subsequently . 

The approximate solutions deviate somewhat from the exact solutions 
for K/fea =» 1 and essentially coincide with them for K/Ra « 1. How- 
ever, near the critical values of Ra -sdiere the approximate solutions 
become infinite, there are, of course, very great discrepancies bet-ween 
■the two solutions. Therefore, solutions -were computed from the exact 
equations for val-ues of Ra at and on either side of the critical. 

These res\ilts are presented in figures 9 and 10 for the case of unequal 
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wall temperatiores (n » 0 ) . It should he recalled that in reference 1 
it was first pointed out that solutions of the exact eguations (taking 
into account the effects of frictional heating) were not unique. In 
fact, two distinct solutions were found for each given set of parame- 
ters up to a limiting set in the only problems which have been solved 
exactly (see refs. 1 and 5). Beyond the certain limiting set of pa- 
rametric values, no exact solutions could be found. The conjecture 
is that this corresponds to a choking condition. 

The two sets of exact solutions for several Bayleigh numbers (see 
table II ) at smd near the first critical value of Ra, shown in figures 
9 and 10 for a = 0 and n = 0, axe denoted by u(^) and u(2). The 
superscripts axe used merely to identify a given exact solution, and 
the (l) is used to denote that solution with the algebraically smaller 
initial slope. Computations from the corresponding approximate solu- 
tions are included in figures 9 and 10. For the case of equal wall 
temperatures (n = l), the velocity extrema of both exact solutions are 
superimposed on figure 6. It can be seen from figures 6, 9, and 10 
that the flows ccorputed from the approximate solutions correspond to 
different exact solutions on either side of the first critical Ra. 

Beyond the first critical Ra. separate consideration must be 
given to the symmetric case (n = 1} and the asymmetric case (n =» O) . 

For n =■ 1, recall that Ra =« «4 ig the only critical Raylei^ num- 
ber (in the range computed, l.e., from Ra =* 10 to 1600) indicated by 
the approximate solutions. Thus, for 97.41 < Ra < 7890.13 (the latter 
being the second critical point for n = l), the approximate so- 

lutions indicate flows of the same general character. However, the 
exact solutions of the superscript (l) type for n = 1 appear to be 
discontinuous at approximately Ra = 1579 (see fig. 6). The change 
in flow type indicated thereby for the larger (than 97.41) Ra may 
iD^Jly that the frictional heating delays (in a Raylel^ number sense) 
the Instability. Rote that the approximate solution always corresponds 
to that exact solution with the lower velocities. 

Exact (see table II ) and approximate solutions near the second 
critical value of Ra ((2rt)^) for a = 0 and n = 0 are shown in 
figures 11 and 12. The two exact solutions are shown only for 
Ra = 1758.55. Although the flow in adjacent sections of the channel 
is in opposite directions, the qualitative effects are essentially 
the same as those near the first critical Ra. No uniisual behavior 
of the exact solutions for n = 0 was foimd, but not as many of these 
solutions for n =« 0 were obtained as for n = 1. 

Near the critical points there are appreciable quantitative and 
qualitative differences between the exact and approximate solutions, 
so ttiat the exact solutions or hl^er-order approximations should be 
used there for reasonable accuracy. Iherefore, it appears that near 
the critical Ra the effects of frictional heating are very ln^ortant. 
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It has "been pointed out (figs. 6 and 9 to 12) that on either side 
of the critical points the solutions neglecting frictional heating cor- 
respond to different exact solutions hut always those with the smaller 
velocity extremum. The question then arises: Since the exact solutions 

do not display a discontinuity (with Ra) at the values indicated hy the 
approximate solutions, is the aotual flow really of a different type on 
either side of the critical Ra as is indicated hy the linearized solu- 
tions? Although a more general study of the present configuration would 
he necessary (perhaps relaxing the condition of fully developed flow) to 
answer this question conclusively, it seems reasonahle that the actual 
flow does in fact change character in accord with the linearized solu- 
tions and the exact solutions for n = 1. The reason for this statement 
is that the approximate solutions are approximate only in that friction- 
al heating is neglected. In starting a flow, the frictional heating ef- 
fects are not of primary importance until velocities of appreciable mag- 
nitude are encountered. Therefore, it is felt that the larger velocity 
flows. Indicated hy one of the two exact solutions, result frcm a regen- 
erative action of the frictional heating (see ref. 5) and that for a 
given set of conditions the flow first established would he that with 
the smaller velocities. Also, the changes in the character of the flow 
on either side of the critical Ra in the experiments with unstable 
horizontal fluid layers seem to lend further support to this contention. 

Near the critical, points themselves, however, the velocities are 
always large; and, therefore, ■ the frictional heating may play an impor- 
tant role in the transition from one type of flow to the other. 

The velocity and temperature distributions for the case of no net 
mass flow in the channel will he affected hy frictional heating in ex- 
actly the same manner as was discxissed above. 


Both Side Walls Insulated 

Ho frictional heatj.ng . - The effect of heat sources in the fluid 
for no frictional heating and with both side wadis insulated can he 
studied hy comparing equation (23) with the solution for a = 0 as 
given in reference 12. Por example, in reference 12 it was shown that 
for no heat sources flow does not ensue except for values of Ra that 

satisfy the equation cos RalAcosh Ral/^ = 1. However, when a 0, 
the present results indicate that the flow and heat transfer are much 
the same as for the linear sxnrflace temperature cases, except that the 
critical values of Ra are hi^er (the first is 3803.22). A repre- 
sentative set of profiles for a = 10, K = 10, and Ra = 10, 7^, and 
8^, with both walls insulated, is shown in figures 13 and 14. Again, 
the different regimes on opposite sides of the critical Ra are 
evident . 
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grlctioEial heating . - Solutions including the effects of frictional 
heating obtained by means of the hi^-speed computing machine for 
Ra = 7^ and 8^ are presented in figin-es 13 and 14. The initial values 
of the exact solutions are presented in table III. These solutions co- 
incide with those computed neglecting frictional heating, as is expected, 
since the ratio K/Ra « 1. At the critical Rayleigh number (3803.22), 
no solution could be obtained with frictional heating. This insolvabll- 
ity may be the result of the "choking condition" previously mentioned 
and discussed in reference 1. 


CONCLUDING REMARKS 

An analysis of the flow subject to a body force between two paral- 
lel plane surfaces oriented parallel to the body force direction with 
heating from below and including the effects of heat sources and fric- 
tional heating shows that the basic physical characteristics are similar 
to flows in horizontal fluid layers with heating fi*om below. The prime 
characteristics in this respect are the existence of critical Rayleigh 
numbers and the difference in the flow and heat transfer on opposite 
sides of critical Rayleigh number. 

Detailed velocity and temperature distributions were computed for 
linearly varying wall temperatures vb.en the wall temperatures were spec- 
ified or when both walls were Insulated, and also for the special case of 
no net mass flow in the channel. The critical Rayleigh numbers are the 
same for corresponding flows with and without heat sources for the spec- 
ified wall tenrperatiAres but are different when both walls are insulated. 
Ccmrparison of two identical configurations, in one of which, however, 
the fluid was heated from below, indicated that the flow and heat trans- 
fer can be appreciably affected by heating from below. 

The effects of frictional heating were once again found to be im- 
portant when K/Ra is of unit order and are also important in the vi- 
cinity of the critical Ra. Two distinct solutions were obtained for 
each set of parametric values when frictional heating is considered. 

The solutions obtained neglecting frictional heating correspond to one 
of the two sets of exact solutions except near the critical regions. 

The approximate solutions on one side of a critical region correspond 
to one of the two distinct solution types obtained with frictional 
heating but on the other side of the critical region correspond to the 
other type exact solution. Hence, the approximate solutions indicate 
completely altered regimes on opposite sides of the critical Rayleigh 
number, analogous to the situation experimentally found for horizontal 
layers heated from below. However, close to the critical region the 
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approximate solutions are invalid; and, hence, the question as to -which 
of the two exact solutions properly describes the flow in this region 
can be ans-wered perhaps by a more general analysis. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, April 29, 1955 
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APPENDIX - SYMBOLS 


The following symbols are used in this report: 
longitudinal ten^jerature gradient, ST*/^X 
parameter defined by eq. (4) 
constants, i = 1, 2, 3, 4 
specific heat at constant pressure 
distance between panadlel plates 

negative of X- component of body force per unit mass 
modified Grashof nimiber, 3fx(-A)d^/v^ 
constant, -Ra^/^CK/64 


Bfyd 

frictional-heating parameter, PrGr 

integers 

dimensionless mass flow 


P^p^f|(-A)d^ 

xp 




wall-temperature parameter, 1 + 


C(-A)d 


presBonre 

Prandtl number, Cpp/x 
heat due to heat sources 
constant, Ra/l6 

modified Rayleigh number, i PrGr 

temperature 

velocity 

d i mensionless velocity 

dimensionless velocity defined by eq. (6) 




3573 
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J 


X 

Y 

y 

a 

3 


longitudinal or axial (vertical) coordinate 
transverse (horizontal) coordinate 
dimensionless transverse coordinate, Y/d 
dimensionless heat-source parameter, Qd/(-A)x 


volumetric expansion coefficient. 



X ratio of specific heats 

T] dimensionless coordinate, 2y - 1 

0 temperature difference, T - or t 

X coefficient of thermal conductivity 

\ constant, oK/64 

li. absolute viscosity coefficient 

V kinematic viscosity coefficient 

p density 

T dimensionless temperatvire difference, KB/(-A)d 

Subscripts ; 

e extremum, value 

■Wq conditions at y = 0 

■w-]_ conditions at y = 1 
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TAHIE I. - EXACT S0LUTI0K8 FOR HALLS AT SAME TEMEE)RATURE 


(a) Ra - 10; a > 0 


[n - 1, K - 10, C - -1.] 

(b) Ra - 10, a » 10 




y 

u(l) 

u(l)’ 


ua)”’ 


y 

u(l) 


u(l)" 

uW" 

0 

0 

- 5.205 


- 0.2362 


0 

0 

- 0.8240 

10.00 

- 50.06 

.10 

-.4704 

- 4.202 

10.09 

1.742 


.10 

-.04033 

-.06760 

5.494 

- 40.07 

.20 

-.8399 

- 3.183 

10.31 

2.452 


.20 

-.02488 

.3081 

1.986 

- 30.10 

.30 

- 1.106 

- 2.140 

10.54 

2.189 


.30 

.01126 

.3729 

- .5240 

- 20.09 

.40 

- 1.267 

- 1.076 

10.72 

1.262 


.40 

.04300 

.2367 

- 2.032 

- 10.05 

.SO 

- 1.521 

.0001010 

10.79 

- .002926 


.50 

.05526 

.00004873 

- 2.534 

- .001519 

.60 

- 1.267 

1.077 

10.72 

- 1.267 


.60 

.04301 

-.2366 

- 2.032 

10.05 

.70 

- 1.106 

2.141 

10.54 

- 2.194 


.70 

.01128 

-.3728 

-.5246 

20.09 

.80 

-.8399 

3.183 

10.30 

- 2,457 


.80 

-.02486 

-.3082 

1,985 

30.10 

.90 

-.4704 

4.202 

10.08 

- 1.748 


.90 

-.04032 

.05746 

5.493 

40.06 

1.00 

- . 1986 X 10 -* 

5.204 

9.997 

.2300 


l.OO 

-. 1445 X 10 '* 

.8237 

9 . 99 S 



50.06 


(c) Ra « 100, a = 0 (d) Ra. » 1600, a » 0 


y 

u(2) 


IQ 

u(2)'" 


7 


u(3-)' 

uW” 


0 

0 

15.34 

10.00 

- 227.2 


0 

0 

- 0.8000 

10.00 

- 31.05 

.10 

1.547 

15.25 

- U .27 

- 195.6 


.10 

-.03526 

.04072 

6.745 

- 34.98 

.20 

2.985 

13.22 

- 28.72 

- 152.0 


.20 


.5321 

3.017 

- 38.72 

.30 

4.140 

9.666 

- 41.48 

- 102.6 


.30 

.05839 

.6431 

-.7247 

- 34.44 

.40 

4.884 

5.090 

- 49.18 

- 51.42 


,40 

.1358 

.4182 

- 3.538 

- 20.33 

.50 

5.141 

. 3334 X 10-3 

- 51.75 



.50 

.1354 

-.004572 

- 4.574 

.1771 

.60 

4.885 

- 5.089 

- 49.18 

51.44 


.60 

.1130 

-.4256 

- 3.505 

20.61 

.70 

4.140 

- 9.665 

- 41.47 

102.7 


.70 

.05699 

-.6461 

-.6719 

34.54 

.80 

2.986 

- 13.22 

- 28.71 

152.0 


.80 


-.5297 

3.068 

58.59 

.90 

1..548 

- 15.25 

- 11.26 

195.6 


.90 

-.03622 ^ 

-.05424 

6.773 

34.65 

1.00 

.6819x10'3 

- 15.34 

10.01 

227.3 


1.00 

-. 2323 > a 0-3 

.8073 

9.988 

50.63 


(a) Ra - 1600, a - 10 


y 

«{ 1 ) 

a(l)' 

u(D” 


0 

0 

iHH 

10.00 

- 39.61 

.10 

-.05429 

-.1666 

6.342 

- 35.00 

.20 

-.04503 

.2944 

2.900 

- 33.56 

.30 

-.006514 

.4236 

-.2216 

- 27.84 

.40 

.03053 


- 2.449 

- 15.71 

.50 

.04432 

- .01711 

- 3.219 

.6815 

.60 

.02733 

-.3076 

- 2.321 

16 . 81 - 

.70 

-.01168 

- .4341 

-.01617 

28.25 

.80 

-.05018 

-.2837 

3.104 

33.11 

.90 

-.05745 ^ 

.1942 

6.464 

33.87 

1.00 

. 2669 X 10 -° 

1.015 

9.991 

38.22 
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TABLE II. - EXACT SOLUnONS FCE UNEQUAL WALL TEMIERATURES 
[n = 0, K = 10, C = -1.] 

(a) Ba = 10, a = 0 (t) Ha = 10, a > 


- .3015 
-.5116 
-.6389 
.40 -.6917 
.50 -.6793 
.60 -.6115 
.70 -.4985 
.80 -.3514 
.90 -.1813 
1.00 .3222JJ10" 


i.555 -10.11 
1.513 -10.71 
5-419 -11.14 
i.292 -11.35 
L.157 -11.33 
.03429 -11.09 


.1571 
.3566 
.5507 
.40 .7021 

.50 .7846 

.60 .7827 

.70 .6916 

.80 .5177 

.90 .2778 

1.00 


1.172 
1.872 
2.039 
1.780 
1.204 
.4195 
-.4660 -8 
-1.344 -8 

-2.107 -6 


.999 

.384 

.682 

.885. 

.009144 


0 0 I 

.10 -.7594 

.20 -1.384 
.30 -1.842 
.40 -2.107 
.50 -2,166 
.60 -2.023 
.70 -1.695 
.80 -1.218 
.90 -.6353 

1.00 .3330¥10~3 


u(l)’ 


-8.150 

20.90 

-6.980 

20.49 

-5.464 

17.64 

-3.647 

12.95 

-1.630 

7.027 

.4405 

,4475 

2.401 

-6.239 

4.092 

-12.49 

5.381 

-17.72 

6.179 

-21.32 

6.446 

-22.65 


u(2)” 


32.42 

35.18 

30.71 

20.23 

5.551 

-11.05 

-27.15 

-40,55 

-49.62 

-53.45 

-51.88 


-297 .1 
-245.2 
-184.3 
-122.9 
-65,17 
-10.67 
43.61 
101.1 
163.1 

227.0 

285.0 


(d) Ea = 97,41, a 


.10 -1.282 
.20 -2.384 
,30 -3.226 
.40 -3.741 
.50 -3'. 886 
.60 -3.655 
-3.077 
-2,216 
.901-1.157 
L.OO .156CP«10'^ 


0 

1.107 

2.160 

3,031 

3.616 

3.850 

3.702 

3.182 

2.332 

1.235 

.291010.0-2 


-6.955 


10.00 

17.97 

26.15 

32.99 

37.20 

37.93 

34.94 
28.59 
19.82 

9.843 I 
-.003217 
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TABLE II. - Continued. EXACT SOLUTIONS FOR UNEQUAL WALL TEMPERATURES 
[n - 0, K =■ 10, C = -1.] 

(e) Ra = 100, a. »• 0 


y 





0 

0 

9.877 

10.00 

-150.7 

.10 

1.013 

10.14 

-4.395 

-135.4 

.20 

1.984 

9.065 

-16,77 

-110.9 

.30 

2.789 

6.883 

-26. .37 

-80.33 

.40 

3.334 

3.899 

-32.73 

-46.42 

.50 

3.554 

.4525 

-35.61 

-11.12 

.60 

3.420 

-3.105 

-34.95 

24.34 

.70 

2.941 

-6.420 

-30.78 

58.81 

.80 

2.156 

-9.149 

-23.27 

90.74 

.90 

1.141 

-10.97 

-12.79 

117.7 

1.00 

.1563«.Ct3 

-11.63 

-.2000x10-3 

136.5 


(f) Ra = 117.41, a - 0 


B 

u ( l ) 

u ( 2 ) 

u ( l )* 

u ( 2 )' 

u(D” 

u ( 2 )” 

„a)'" 

u ( 2 )'" 

0 

0 

0 

- 22.65 

5.400 

10.00 

10.00 

129.2 

- 98.41 

.10 

- 2.192 

.5738 

- 20.93 

5.916 

24.95 

.4285 

164.3 

- 91.80 

.20 

- 4.132 

1.153 

- 17.60 

5.520 

41.69 

- 8.128 

164.6 

- 78.27 

.30 

- 5.657 

1.652 

- 12.65 

4.345 

56.70 

- 15.04 

130.1 

- 59.19 

.40 

- 6.620 

2.002 

- 6.427 

2.582 

66.74 

- 19.84 

66.97 

- 36.31 

.50 

- 6.921 

2.156 

.4536 

.4578 

69.56 

- 22.24 

- 11.89 

- 11.42 

.60 

- 6.533 

2.090 

7.217 

- 1.781 

64.40 

- 22.11 

- 89.65 

13.81 

.,70 

- 5.507 

1.805 

13.10 

- 3.882 

52.21 

- 19.52 

- 150.1 

37.73 

.80 

- 3.963 

1.326 

17.50 

- 5.608 

55.34 

- 14.67 

- 182.1 

58.59 

.90 

- 2.067 

.7023 

20.11 

- 6.753 

16.88 

- 7.962 

- 181.7 

74.51 

1.00 

- . 4120x10-3 

. 8901XL0 -^ 

20.93 

- 7.159 

- 002140 

.007880 

- 151.5 

83.61 


(g) Ra = 1000, a = 0 


(h) Ra = 1000, a - 10 


B 

u(l) 

u(l)’ 

u(l)" 


0 

0 

-2,040 

10.00 

8.015 

.10 

-.1528 

-1,007 

10.52 

-.2448 

.20 

-.2017 

.01490 

9.608 

-18.79 

.30 

-.1561 

.8490 

6.762 

-37.35 

.40 

-.04417 

1.317 

2.419 

-47.62 

.50 

.09160- 

1.319 

-2.330 

-45.07 

.60 

.2048 

.8815 

-6.165 

-29.75 

.70 

.2581 

.1536 

-7,995 

-5.961 

.80 

.2336 

-.6327 

-7.308 

19.29 

.90 

.1379 

-1.231 

-4.340 

38.47 

1.00 

-.1136X10-^ 

-1.455 

-.01511 

45.74 


D 

u ( l ) 

u ( l )' 

u ( l )” 

u ( l )”’ 

0 

0 

- 2.412 

laoo 

2.010 

.10 

-.1906 

- 1.394 

10,36 

2.002 

.20 

-.2783 

-.3673 

9.927 

- 12.22 

.30 

-.2682 

.5339 

7.793 

- 30.29 

.40 

-.1815 

1.137 

4.045 

- 43.20 

.50 

-.05509 

1.317 

-.4729 

- 45.02 

.60 

.06701 

1.057 

- 4329 ' 

- 34,06 

.70 

.1452 

.4663 

- 6.941 

- 12.89 

.80 

.1560 

-.2500 

- 6.955 

12.78 

.90 

.09944 

-.8405 

^.465 

36.09 

1.00 

-. 1579X10 “^ 

- 1.077 

-.009459 

51.35 
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TABLE IT. - Concluded. EXACT SOLUrrONS FOR UNEQUAL HAIX TEltPERATURES 
[n = 0, K = 10, C = -1.] 

(i) Ra = 1358.55, a - 0 (j) Ha = 1558.55, a 


0 0 
.10 -.3982 
.20 -.6003 
.30 -.5592 
.40 -.3035 
.50 .06772 

.60 .4174 

.70 .6188 

.00 .5991 

.90 .3650 

1.00 -.7953X10" 


-4.660 

-3.135 

-.8246 

1.598 

3.346 

3.841 

2.934 

.9675 

-1.349 

-3.188 

-3.887 


u(l) 

tit 

109, 

.4 

82. 

.18 

12. 

.18 

-69. 

.28 

-129. 

.2 

-144. 

.4 

-109. 

.2 

-36. 

.11 

49. 

.30 

117. 

.5 

144. 

.4 


) 0 
.10 -6 .005 
.20 -9.854 
.30 -10.15 
.40 -6.837 

.50 -1.309 
.60 4.262 

.70 7.809 

.80 8.096 

.90 5.094 


-64.10 

-52.45 1 221.0 
-22.061 367.9 
16.32 


•6.837 47.44 
•1.309 59.33 
4.262 48.52 
7.809 20.29 
8.096 -14.59 
5.094 -43.41 
-.1463x10-5 -54.84 


-755 .5 
-2002 
-2345 
-1789 
-679.4 
616.7 
1779 
2452 


1600, a 


(1) Ra - 1600, a - 10 


0 

1.800 
2.938 
2.944 
.40 1.795 
.50 -.08615 
.60 -1.982 
.70 -3.155 
.80 -3.140 
.90 -1.934 


18.88 

15.79 

6.163 

-6.076 

-16.18 

- 20.20 

-18.44 

- 6.210 

6.530 

16.78 


1.00 -.7300x10-5 20.64 


10 

68 

17 


-119 .5 
-65.74 
-1.733 
74.70 

123.0 

123.1 
75.13 


-8.54.1 

-673.5 

-267.8 

220.5 

627.3 
804.9 
670. .5 

260.4 
-259.0 
-663.2 


-.01997-785.5 


y 

u(2) 

u(2)’ 

u(2)" 

u(2)"' 

0 

0 

18.81 

10.00 

-866.8 

.10 

1.791 

15.68 

-69.51 

-677.1 

.20 

2.913 

5.959 

-117.9 

-264.3 

.30 

2.896 

-6.334 

-119.7 

228.2 

.40 

1.721 

-16.42 

-75.12 

635.7 

.50 

-.1793 

-20.34 

- .3720 

810.6 

.60 

-2.082 

-16.42 

76.36 

670.8 

.70 

-3.244 

-6.030 

124.4 

255.2 

.80 

-3.204 

6.821 

123.8 

-266.5 

.90 

-1.967 

17.11 

75.23 

-667.8 

1.00 

-.1390X10"^ 

20.97' 

-.934ttcl0"3 

-781.5 


m) Ra = 1758.55, a 


0 

.3381 
.5909 
.6309 
.4286 
■ .05856 
-.3334 
-.5902 
-.6084 
-.3809 
-.6685X10-4 


0 

14.95 
22.32 
18.05 
2.146 
-21.13 
-43.32 
-54 .42 
-48.51 
-27.78 

.2300X10-4 


uU)' u(2)’ 


163.9 
121.7 

18.91 
-104.3 
-206.8 
-244 .3 
-181.8 
-29.47 
143.5 

256.9 
287.0 


u(2)" B u(2)"’ 


10.00 
-9.026 
-22.14 
-25.15 
-17.41 
-2.2R6 
15.99 
24.71 
25.45 
15.87 
-. 06536 


10.00 

-792.1 

-1196 

-1200 

-775.1 

91.46 

1145 

1777 

1535 

690.6 


-200.7 

-169.8 

-85.07 

26.06 

122.9 

168.7 

145.1 

61.66 

-47.51 

-137.1 


-.005290 -170.4 -3740 
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TABLE III. - INITIAL VALUES OF EXACT SOLUTIONS FOR BOTH WALLS INSULATED 


[oK = 0.01.] 

(a) Ra = 2401 


u^l^(O) 



u(2)’(o) 

yW"(0) 




0 

0 

-0.000035 

166.0 

-0.0004495 

-3172 

0 

0 


(B) Ra = 4096 



u(^^{0) 

y(^^’(0) 

u(2)'(0) 


u(2)”(o) 



0 

0 

-50 

-0.00028 

387.1 

-1 

0.002085 

0 

0 
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Figure 1. - Schematic sketch of configuration considered. 
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< 0 

-1 




o1 

lOJ 


> 0 

1 


ms 




(a) Rayleigh number, 10. 


(b) Rayleigh number, 61 


(o) Rayleigh number, IQO. 


(d) Rayleigh number, 1600. 


Figure 2. - Velocity dletributlona with and without heat souroea Tor heating from below 
Ta < 0) and not (A> 0). Identical wall temperaturee (n - l)j K - 10. 
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(o) Rayleigh xxunbery 


(d) Rayleigh mimbery 1600. 


Figure 5. - Teb^erature dlstributlona with and without heat sources for 
heating from below (A< 0) and not, (A > O). Identical wall temperatures 
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(a) Rayleigh nunber, 10. 


(b) Rayleigh number, 100. 



(c) Rayleigh number, 1000. (a) Rayleigh number, ISOO. 

Figure 5. - Temperature dlatribution with and without heat aources for heating 
from below (A < 0) and not (A > O). Bnequal wall temperaturea (n - 0); K - 10. 
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Figure 9. - Velocity distributions for Rayleigh numbers at and on 
both sides of first critical Raylei^ number. Unequal wall tem- 
peratures (n = O); K = lOj a ™ Oj C = -1. 
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y = T/d 


Figure 10. - Teimperature dlstrlTjutions for Rayleigh numbers at and 
on both sides of first critical Rayleigh number. Unequal wall 
temperatures (n »• 0)j K = lOj a = Oj C = -1. 
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y - Y/d 

Figure 11. - Velocity distributions for Kaylel^ numbers at and on both sides of 
second critical Rayleigh number. Uneiual wall temperatures (n - O); K = 10} a = 0; 
C - - 1. 
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Figure 12. - Temperature distrilsutions for Rayleigli numl^ers 
at and on "both sides of second critical Rayleigh number. 
Unequal vail temperatures (n=0)} K= 10} a=0) C = -1. 







